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ABSTRACT 

In this paper, predictions are made of noise radiation 
from single, supersonic, axisymmetric jets. We exam- 
ine the effects of changes in operating conditions and 
the effects of simulated enhanced mixing that would in- 
crease the spreading rate of the jet shear layer on radi- 
ated noise levels. The radiated noise in the downstream 
direction is dominated by mixing noise and it is well 
described by the instability wave noise radiation analy- 
sis. A numerical prediction scheme is used for the mean 
flow providing an efficient method to obtain the mean 
flow development for various operating conditions and to 
simulate the enhanced mixing. Using far field radiated 
noise measurements as a reference, the calculations pre- 
dict that enhanced jet spreading results in a reduction 
of radiated noise. 

1 INTRODUCTION 

Recent interest in the development of a new High 
Speed Civil Transport has renewed research efforts to 
make the aircraft environmentally compatible. One is- 
sue is the level of noise generated by the supersonic jets 
exhausting from the propulsion systems. Supersonic jets 
axe intense noise generators and means must be found to 
modify the noise generation process to reduce radiated 
noise levels in order for the aircraft to meet commu- 
nity noise regulations. A recent review by Seiner and 
Krejsa [1989] discusses the challenge of reducing super- 
sonic jet noise associated with both mixing and shocks 
while maintaining acceptable propulsion system perfor- 
mance requirements. In this paper, we focus on the issue 
of mixing noise and the concept that the promotion of 
rapid mixing of the jet shear layer is an effective means 


of reducing radiated levels of mixing noise. 

At present, the primary means of promoting enhanced 
jet shear layer spreading is either through acoustical or 
mechanical excitation. It is well known that subsonic 
jet mixing can easily be enhanced through acoustic ex- 
citation; but, for supersonic jets, acoustic excitation is 
much less effective (Lepicovsky et al. [1985]). Mechani- 
cal means for promoting jet shear layer spreading were 
found to be more effective in supersonic jets than acous- 
tic excitation. For example, Ahuja and Brown [1989] 
inserted small tabs into the jet flow at the lip of the 
nozzle exit. Alternatively, Gutmark et al. [1989] used 
nozzle shaping, such as elliptic and rectangular nozzles, 
to promote increased shear layer spreading. Experi- 
ments to measure the noise from supersonic jets with 
mechanical tabs have been performed by Ahuja [1993] 
and Kobayashi et al. [1993]. In both studies, the jets 
contained shocks and for the most part the overall noise 
reduction was primarily due to the reduction of noise as- 
sociated with shocks. However, the studies did show de- 
creases in noise levels in the downstream direction. Mea- 
surements of noise from perfectly expanded jets from el- 
liptic nozzles (Seiner and Ponton [1991]) showed noise 
reductions compared to the equivalent axisymmetric jet. 
In this case, the results depended on the azimuthal an- 
gle. 

While there have been noise predictions for both su- 
personic axisymmetric and elliptic jets, noise predictions 
have not been made for supersonic jets with enhanced 
mixing. It is the purpose of this paper to make noise 
predictions for single, supersonic, axisymmetric jets for 
given operating conditions and evaluate the resulting 
noise reduction from simulated enhanced mixing that 
increases the spreading rate of the jet shear layer. 
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It is known that instability waves are the dominant 
source of mixing noise radiating into the downstream 
arc of a supersonic jet when the phase velocities of the 
wave are supersonic relative to ambient conditions. The 
analysis to predict the noise radiated to the far field is 
based on the determination of the axial growth and de- 
cay of the instability wave in the jet shear layer. In order 
to complete the analysis, the mean flow properties are 
needed. In the next section, we describe the numerical 
prediction method for the mean flow development. This 
provides am efficient method to obtain the mean flow 
for various operating conditions and to simulate shear 
layer growth due to enhanced mixing. The procedures 
to calculate the evolution of instability waves and their 
radiated noise are presented next. This is followed by 
a comparison of calculated far field radiated noise pat- 
terns to measured supersonic jet data in order to cali- 
brate one unknown constant. Based on those results, we 
then compare calculated results for a jet with simulated 
enhanced mixing to measured data. Finally, we present 
a numerical study of the changes in radiated noise levels 
with changes in jet spreading rate. 

2 MEAN FLOW PREDICTIONS 

This section discusses the numerical procedure used to 
calculate the mean flow development of a compressible, 
axisymmetricjet. The jet flow is assumed to be perfectly 
expanded, thus the jet static pressure is matched to the 
ambient pressure. Often in the analysis of the mixing 
noise generation in single supersonic axisymmetric jets, 
the axial variation in the mean velocity profile is based 
on simple analytic functions with axially varying param- 
eters described by equations fitted to measured data. In 
supersonic jets where temperature profiles are not easily 
measured, Crocco’s relation is used to obtain mean den- 
sity profiles. However, in order to predict radiated jet 
noise for a variety of operating conditions, it is appro- 
priate to use numerically generated mean flow profiles. 
Thus, we use a set of compressible, Reynolds averaged, 
boundary layer equations with a modified mixing length 
model to determine the Reynolds stresses as the basis 
for the numerical mean flow analysis. To simulate the 
effects of enhanced mixing, the coefficient of the mixing 
length model is increased to provide a larger eddy vis- 
cosity which causes the jet shear layer to spread faster. 

2.1 Turbulence Model 

The development of a prediction scheme for the mean 
flow properties of a jet was completed with the capabil- 
ity to compute both single stream axisymmetric jets and 
dual stream coaxial jets, Dahl [1994]. Simplicity and ro- 
bustness were emphasized in order to calculate the mean 
velocity and density and their derivatives accurately for 
later usage in the instability wave calculations. This led 
us to choose a simple turbulence model, resulting in a 
high level of empiricism. 

The compressible equations of motion are simplified to 


the boundary layer form. The assumption is also made 
that the density- velocity correlations may be neglected. 
For the single jet case, the Reynolds stress and the heat 
flux terms are described by the following mixing-length 

model: _ 
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l is a characteristic mixing length scale given by 
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and k is a constant to simulate the effects of enhanced 
mixing. 

#c = 1, no enhanced mixing 

k > 1, enhanced mixing ^ ' 

The factor C\ is the incompressible part of the mixing 
length constant. It depends on the the velocity ratio 
Uco/Uj and the density ratio Poo/p.>* The Ci factor is 
the compressible part of the mixing length constant. Its 
purpose is to decrease the growth of the shear layer as 
compressibility effects become important. It depends 
on a Mach number in a frame of reference convecting 
with the real phase speed of a growing disturbance in 
the shear layer. This convected Mach number depends 
on the velocity ratio, the density ratio, and the Mach 
number of the jet. Thus, both factors depend on the flow 
conditions and a calibration procedure was developed to 
obtain empirical equations to describe both factors. The 
details are given in Dahl [1994]. 


2.2 Numerical Method 


The numerical method initially follows the stream 
function approach given by Crawford and Kays [1976]. 
The equations of motion in boundary layer form are 
transformed into stream function coordinates using 


A ~ m f*\ 
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where v is the mass averaged radial component of the 
mean velocity. The boundary layer equations become 
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where H is the total mean enthalpy and 


= P + PT 
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The numerical equations are derived by using fully im- 
plicit differencing on equations (7) and (8). This insures 
that the numerical problem is inherently stable. Since 
both the mean flow and instability wave problems are 
calculated using the same grid, maintaining sufficient 
grid resolution required that the problems be solved on 
a fine, evenly spaced r-grid. Details of the finite differ- 
ence algorithm and solution procedure are given in Dahl 
[1994] as are favorable comparisons between calculated 
and measured data for various subsonic and supersonic 
jets. 

3 INSTABILITY WAVES AND RADIATED NOISE 

It is well known that thin free shear layers containing 
an inflection point in the mean velocity profile are inher- 
ently unstable even in the absence of viscosity. An insta- 
bility wave in the shear layer initially grows rapidly. As 
the shear layer grows, the wave growth rate decreases. 
Eventually, the shear layer is too thick to support un- 
stable waves and the wave amplitude decreases until it 
disappears. This instability wave process is assumed 
to be governed by the linearized, inviscid, compressible 
equations of motion. 

For slowly diverging jet flows, two solutions are cre- 
ated that apply to separate but overlapping regions. Fol- 
lowing the approach of Tam and Burton [1984], the inner 
region, including the jet flow and the immediate region 
just outside the jet, has different length scales between 
the radial and the axial directions that leads to a mul- 
tiple scales expansion of the governing equations. With 
the pressure disturbances represented as 


the slow mean flow development in the axial direction 
s = ex , ^(s) is an axial phase function related to the ax- 
ial wavenumber or by d<f>/ds = a($), n is the azimuthal 
mode number, and u> is the radian frequency. 

The general solution to equation (12) is written 

po(r, s) = A 0 (s)Cf(r, s) + Bo(*)C 3 P (r, s). (13) 

As r — ► oo, the flow outside the jet becomes uniform 
and the solution to equation (12) may be written in the 
form 

Po = A 0 H^(iX(a)r) + B 0 H™(i\{a)r) (14) 

where 

A(a) = [a 2 - - aUoo) 2 ] 1/2 (15) 

and H ^ and H ^ are nth-order Hankel functions of 
the first and second kind, respectively. Thus, for large 
r, equation (13) must tend to equation (14). As r — ► 0, 
equation (13) must be finite. 

In the outer region, which slightly overlaps the in- 
ner region, the governing equations control disturbances 
that are acoustic in nature. These disturbances have the 
same length scales in all directions; hence, all coordi- 
nates are treated equally. To create am outer solution in 
a form that allows it to be asymptotically matched to the 
inner solution, we use the axial coordinate s = ex and 
the radial coordinate r = er. The solution is obtained by 
Fourier transforming the outer region governing equa- 
tions in the s direction. After considerable algebra, the 
lowest order outer solution is 

oo 
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the lowest order set of equations may be reduced to 
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which is known as the compressible Rayleigh Equa- 
tion. In equation (11), £ m (e) are the gauge functions 
of the aymptotic expansion in the small parameter e 
where 6o(e) = 1, s is the slow variable to recognize 


The two solutions are matched asymptotically in an 
overlap region. To lowest order, we find that 

A 0 (s) = A 0 (s) and £o(*) = 0 . (18) 

With the condition of finite value at r = 0, equation 
(12) has become an eigenvalue problem with solutions 
only for certain values of the eigenvalue a. We have 
used a finite-difference approximation to discretize the 
problem. The eigenvalue is found from the resulting 
diagonal matrix using a Newton-Raphson iteration for 
refinement. 

Since the rate of spread of the jet is slow for high speed 
jets and z is very small, Ao(zx) in equation (17) is taken 
to be constant. Furthermore with c* (a;) found from the 
solution of equation (12) at every axial location, the 


3 



axial phase function is found from <f>{ex)fe = J a(x) dx. 

o 

We can then solve for g(i}) in equation (17), the Fourier 
transform of the instability wave, and subsequently the 
near field pressure disturbance is found from equation 
(16). To obtain the pressure in the far field, equation 
(16), in spherical coordinates, is approximated by the 
method of stationary phase. The resulting sound power 
radiated per unit solid angle is 


DW = $\p\ 2 R 2 = 2 


IgO ?)! 2 

[l - Af*, sin' 




(19) 


The stationary point is given by 

pi£ 2 MjU) cos PooM? «ooW 



and ^ is the polar angle. 

To validate the mean flow and noise radiation predic- 
tion schemes, Dahl [1994] compared calculated results 
using numerically generated mean flow profiles to mea- 
sured data (Seiner and Ponton [1985]) and to calculated 
results using analytically generated mean flow profiles 
(Seiner et al. [1993]) and found favorable comparisons 
in the far field radiated noise patterns even though the 
numerically generated mean flows had longer potential 
core lengths than the measured data. However, the cal- 
culated initial spreading of the jet was similar to mea- 
sured data and encompassed the region of maximum 
growth of the instability wave. Thus, we assume that 
given the jet operating conditions, we can calculate the 
initial jet spreading with sufficient accuracy to obtain 
reasonable results from the instability wave noise radia- 
tion analysis. 

4 PREDICTED AND MEASURED SOUND 

4.1 Variations in Operating Conditions 

Figure 1 shows a comparison between calculated and 
measured far field noise radiation patterns for three 
Mj = 2, hot jet cases given in Seiner et al. [1992]. The 
calculated results are for the first helical mode instabil- 
ity wave at Strouhal number 0.11. The calculations are 
normalized to the measured peak level. By doing this, 
it is found that the constant in the problem, the initial 
amplitude of the instability wave Aq in equation (17), 
increases with jet temperature. This agrees with the 
calculated results of Tam and Chen [1993]. The figure 
indicates that jet spreading rate is changing with the 
operating conditions; but, that far field radiated noise 
at this Strouhal number is increasing with jet tempera- 
ture. The results show that to increase spreading rate by 
changing the operating conditions does not necessarily 
lead to a decrease in radiated noise. As the jet temper- 
ature increases, the convected Mach number increases 
resulting in a shift of the measured noise radiation peak 
to larger angles. This trend is followed in the calculated 


Figure 1: Comparison of Mj = 2 jet calculations to 
measured far field noise data, n = 1, St = .11 


results. 

4.2 Effects of Enhanced Mixing 

One method to enhance mixing is to place tabs around 
the nozzle lip. Kobayashi et al. [1993] ran a series of 
tests using various size tabs. Though their jets are un- 
derexpanded, the noise in the downstream direction is 
dominated by mixing noise and is affected by changes 
in jet spreading. In the experiment, two small tabs with 
0.5% blockage are placed at the nozzle lip. Even though 
the tabs are small and they eliminate jet screech, the 
two tab arrangement could still distort the jet flow ax- 
isymmetry (See Samimy et al. [1993]). Hence, we use 
the measured noise data as a reference in our enhanced 
mixing simulation. The measured OASPL directivity 
patterns for a jet with and without tabs are shown in 
Figure 2. 

To simulate enhanced mixing, k in equation (3) is in- 
creased above 1. Assuming that the initial amplitude 
of the instability wave remains unchanged, the level of 
enhancement is increased until the enhanced jet far field 
noise radiation pattern approximately matches the mea- 
sured data. Figure 2 shows the calculated radiation pat- 
terns for normal and enhanced spreading rate jets at 
Strouhal number 0.20 where the instability wave is near 
its largest total growth. The calculated results are nor- 
malized to the measured data at 20 degrees. Agreement 
was found for k = 1.05 resulting in a 10% increase in 
initial jet spreading rate db/dx where b is the jet shear 
layer half-width. The measured data show a reduction 
in peak level of 4 dB and the calculated data show a re- 
duction in peak level of 3 dB at Strouhal number 0.20. 
In both cases, the peak direction remained unchanged. 
The enhanced spreading rate of the jet decreases the 
maximum amplitude of the growing instability wave re- 
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Figure 2: Effects of enhanced mixing on radiated noise, 
Mj = 1.484. Data: OASPL, Calc.: n = 1, St = .20 


suiting in lower levels of radiated noise. 

Another means of enhancing jet mixing is to use noz- 
zle shaping such as the elliptic nozzle referred to earlier. 
Noise radiation predictions for elliptic nozzles are be- 
yond the scope of this paper, but we use measurements 
of far field noise radiated from elliptic jets as a basis for 
determining how much enhanced mixing an axisymmet- 
ric jet requires to achieve similar noise levels. Again, the 
measured data is in terms of OASPL. We make the cal- 
culations for a Strouhal number of 0.40. Figure 3 shows 
the measured far field noise radiation pattern for the ref- 
erence axisymmetric jet and the elliptic jet on both the 
major and minor axes. The calculated data are refer- 
enced to the peak level of the measured axisymmetric jet 
data. The enhanced jet with k = 1.05 has reduced noise 
levels comparable to the elliptic jet along the minor axis, 
about 3 dB. This is the same level of enhanced spread- 
ing and noise reduction achieved with tabs in Figure 2. 
To achieve comparable levels along the major axis, the 
axisymmetric jet must be enhanced with about k = 1.10 
or about a 21% increase in the initial jet spreading rate. 
The peak noise reduction is about 5.5 dB. 

We next compare the effects of changes in spreading 
due to changes in operating conditions versus enhanced 
spreading. Figure 4 shows the measured and calcu- 
lated data for the two hottest jets shown in Figure 1, 
ie. TJTa = 4.004 and T 0 /T a = 5.422. Even though 
the initial spreading rate has increased from the cooler 
jet to the hotter jet, the effects of increased tempera- 
ture, increased velocity, and increased convected Mach 
number have been to change directivity and slightly in- 
crease the radiated noise levels. If instead of changing 
operating conditions, the spreading rate of the cooler 
jet is enhanced to the same spreading rate as the hotter 
jet, the radiated noise levels are slightly reduced and the 



Figure 3: Effects of enhanced mixing on radiated noise, 
Mj = 1.5. Data: OASPL, Calc.: n = 1, St = .40 

directivity unchanged. This result again assumes that 
the initial instability wave amplitude remains unchanged 
during enhancement. 

Finally, we consider the effects of simulated enhanced 
mixing as a function of spreading rate and Strouhal 
number for a hot, Mj = 2 jet. Figure 5 shows changes in 
peak far field radiated noise levels for enhanced mixing 
relative to the normal spreading rate jet. Enhancing the 
jet spreading rate is more effective at lowering radiated 
noise at lower Strouhal numbers than higher Strouhal 
numbers. Also, the relative increases in spreading rate 
are less effective at noise reduction as the rate increases. 

At present, it is difficult to enhance jet mixing without 
significantly altering the flow making the present analy- 
sis technique difficult to apply. The 10 to 20% increases 
in jet spreading rate we simulated are assumed to take 
place with minimal impact on mean flow profiles. Thus, 
our analysis could proceed showing calculated results 
that are comparable to reference measured data. For 
larger amounts of induced mixing and larger jet spread- 
ing rates, Kobayashi et al. [1993] showed no increases 
in noise reduction. Further studies are necessary to de- 
termine the limits of noise reduction achievable through 
enhanced mixing. 

5 SUMMARY 

In this paper, we have described a procedure for the 
prediction of noise from the instability waves of super- 
sonic, axisymmetric jets. A numerical method was used 
to generate the developing mean flow. This provided 
an efficient means for obtaining the developing mean 
flow for various operating conditions and to simulate en- 
hanced jet mixing. Enhancing the jet spreading rate at 
constant operating conditions decreased the maximum 
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Figure 4: Effects of operating conditions versus en- 
hanced mixing on radiated noise, Afy =2. n = 1, 
St = .11 

amplitude of the growing instability wave resulting in 
decreased far field radiated noise. A change in operat- 
ing conditions to increase the spreading resulted in noise 
increases compared to noise decreases for enhanced mix- 
ing to the same spreading rate. 
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